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The  goals  of  this  project  were  to  attempt  to  simulate  a  fully  viscous, 
three-dimensional  time  accurate  hypersonic  inlet/engine  combined 
geometry,  for  the  specific  purposes  of  studying  the  fundamental  physics 
associated  with  the  phenomena  of  engine  unstarts.  Particular  interest  was 
initially  placed  on  the  interaction  between  the  inlet  shock  structure  and  the 
viscous  boundary  layer,  especially  in  examining  the  growth  of  separation 
regions  and  the  communication  of  downstream  perturbations  to  the 
upstream  inlet  shock  structure. 


ActiviUes  for  the  first  year  were  primarily  focused  on  laying  the 
groundwork  for  accomplishing  the  ultimate  goals  of  this  investigation. 
Required  hardware  and  software  was  obtained  and  implemented.  Grid 
generators  were  tested  and  one  was  selected,  as  was  the  primary 
computational  tool.  On  the  technical  side,  a  substantial  literature  search  was 
conducted.  Several  papers  with  hypersonic  inlet  and  engine  flowfislds  were 
located,  and  have  been  used  for  code  validation.  A  baseline  geometry  was 
selected  for  both  2-D  and  3-D  calculations,  based  on  the  availability  of  data  for 
comparison. 


Gridding  was  completed  for  these  geometries.  A  detailed  study  of  the 
time  accuracy  of  the  code,  vital  to  this  project,  was  undertaken.  Two 
validation  cases  were  considered:  a  quasi-one-dimensional  shock  tube,  for 
validating  internal  flows  (i.e.  combustor  flowfields)  and  a  24°  compression 
corner,  for  validating  external  flows  (i.e.  forebody  inlets).  Extensive 
comparisons  between  the  GASP  results,  analytical  solutions,  and  other  time- 
accurate  codes  were  made. 
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Finally,  time  accurate  studies  on  the  flow  inside  the  combustor  channel 
were  begun  in  the  first  period.  Preliminary  work  examined  the  separation 
region  associated  with  an  unexpected  (off-design)  shock-boundary  layer 
impingement  inside  the  engine  channel,  and  the  resulting  upstream 
propagation  of  that  disturbance. 


Activities  in  the  scond  year  were  directed  towards  beginning  the  three- 
dimensional  Euler,  then  Navier-Stokes  calculations,  to  establish  a  steady-state 
time  accurate  baseline,  which  could  then  be  perturbed  to  study  the  influence 
of  downstream  disturbances.  Both  Euler  and  laminar  Navier-Stokes 
solutions  were  calculated,  though  problems  with  an  ever-expanding  ~ 
separation  zone  in  the  laminar-flow  regime  prevented  the  establishment  of — 
steady-state  conditions. 


I 
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Hardware  and  Software  Support 


In  support  of  the  technical  accomplishments  of  this  work,  the 
following  tasks  were  completed; 
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•  Purchasing  (with  university  funds)  and  installing  a  Silicon  Graphics 
Indigo  Elan  workstation  which  is  dedicated  for  this  project.  Substantial  time 
was  invested  in  setting  up  that  system  and  obtaining  and  installing 
appropriate  software. 

•  Obtaining  several  grid  generation  packages,  including  EAGLE  and 
EAGLEVIEW  from  Mississippi  State  University,  and  becoming  familiar  with 
the  operation  of  these  codes. 


•  Installing  and  becoming  familiar  with  Walter's  GASP  code  as  the 
primary  computational  tool,  and  modifying  that  code  to  facilitate  time- 
accurate  solutions. 

•  Installing  and  becoming  familiar  with  a  graphics  package  (the  FAST 
code  from  NASA  Ames)  for  data  presentation. 


A  major  portion  of  this  research  project,  run  in  parallel  with  the 
validation  and  numerical  experimentation  with  the  computational  tool,  was 
the  selection,  implementation,  and  familiarization  with  a  three-dimensional 
grid  generator.  To  this  end,  several  grid  generation  programs  were 
considered,  including  GRIDGEN,  EAGLE,  and  EAGLEVIEW  (an  interactive 
EAGLE  derivative.) 

The  batch  version  of  the  EAGLE  code  was  first  used  to  generate  a  typical 
3-D  computational  grid  in  the  Langley  3-D  module  geometry.  Figure  21  is  a 
representation  of  one  such  result.  The  wall  geometry  itself,  extracted  from 
the  grid,  is  shown  in  Figure  22.  Note  the  sidewall  compression,  and  the  fact 
that  a  foreplate  is  used  to  produce  an  inlet  boundary  layer.  For  the  work  of 
interest  here,  final  results  will  likely  be  performed  with  the  inlet  forebody 
canted  at  a  more  representative  angle,  or  with  a  preset  boundary  layer  profile. 
One  of  the  ultimate  interests  in  this  study  is  the  inclusion  of  boundary  layer 
effects,  and  the  coupling  of  these  to  the  unsteady  shock  behavior,  so  inclusion 
of  a  realistic  forebody  flowfield  is  of  importance.  Euler  flow  calculations  were 
performed  in  a  two  dimensional  cross-section  of  the  above  geometry  for  a 
baseline  indication  of  the  associated  shock  structure.  Figure  23  shows  density 
contours  and  variations  along  the  centerline  through  the  engine. 

The  time  response  of  the  2-D  RPI  combustor  flowfield  was  calculated  with 
GASP,  using  a  251  x  101  grid  covering  the  inlet,  combustor,  and  nozzle  region 
of  that  engine  shape,  shown  in  Figure  24.  That  inlet  is  not  designed  with  a 
shock-on-lip  configuration,  so  a  strong  shock  reflects  into  the  engine  channel 
and  strikes  the  combustor  wall  boundary  layer.  This  interaction  generates  a 
separated  region  which  propagates  upstream  from  the  intersection  point. 

Figure  25  shows  the  density  and  Mach  number  contours  in  this  engine 
at  freestream  Mach  number  6.  The  relatively  strong  reflected  shock  is  visible, 
as  is  the  thick  boundary  layer.  For  these  calculations  the  boundary  layer  was 
kept  laminar,  though  in  reality  it  will  likely  be  turbulent  past  the  ramp;  a 
turbulent  boundary  layer  will  exhibit  a  much  smaller  separation  region  with 
a  shock  interaction,  so  the  separated  region  is  exaggerated. 


Code  Selection  and  Validation 


For  purposes  of  studying  the  fundamental  physics  of  this  problem,  it  was 
necessary  to  locate  a  computational  tool  which  wap  time-accurate,  three- 
dimensional,  and  had  features  for  including  chemical  models.  To  this  end, 
several  existing  codes  were  examined.  The  TVD  scheme  of  Yee,  et.  al.  has 
been  used  successfully  by  the  principal  investigator  for  both  2-D  and 
axisymmetric  calculations  but  extension  to  three-dimensions  would  have 
represented  a  monumental  effort.  Similarly,  CFL3D  has  also  been  used  at 
Maryland,  but  does  not  include  chemistry,  and  so  was  not  adopted. 

The  General  Aerodynamic  Simulation  Program  (GASP)  of  Walters,  et.  al. 
was  chosen  because  it  had  all  of  the  required  features,  and  was  advertised 
(though  not  validated)  as  time-accurate.  The  first  task  in  using  GASP,  after 
becoming  familiar  with  its  operation,  was,  therefore,  to  validate  the  time- 
accuracy,  after  modifying  the  code  to  facilitate  restart  and  periodic  output. 

Two  cases  were  considered: 

•  Internal  hypersonic  flow  in  a  shock  tube  (one-dimensional  flow,  but 
calculated  in  two  dimensions.) 

•  External  hypersonic  flow  on  Holden’s  24-degree  compression  corner. 

The  shock  tube  calculation  is  shown  in  Figure  1,  which  presents 
density  variation  along  the  tube  as  a  function  of  distance  along  the  centerline, 
for  a  given  pressure  ratio,  using  Roe's  flux  difference  splitting.  This  is 
compared  to  the  exact  Riemann  solution  at  .5  ms.  Note  the  excellent 
agreement  between  the  results.  Temperature  and  velocity  are  shown  in 
Figures  2  and  3,  and  they  also  compare  quite  well  with  the  exact  solutions. 
Corresponding  thermodynamic  contours  are  shown  in  Figure  4 

Subsequent  studies  were  done  examining  different  time-integration 
schemes.  Figure  5  shows  a  comparison  of  two-factor  approximate 
factorization  and  a  Jameson  m -stage  Runge-Kutta  method.  The  implicit 
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calculation  was  run  at  a  very  low  CFL  number  (0.20),  and  agreement  is 
excellent.  The  time  step  was  also  varied  to  establish  the  solution 
independence  over  equivalent  periods  of  time.  Figure  6  presents  the  density 
profile  at  the  same  elapsed  time  (.5  ms)  using  three  different  time  steps,  and 
agreement  is  excellent.  Studies  were  also  done  comparing  Van  Leer  Flux 
Vector  Splitting  and  Roe  Flux  Difference  splitting  algorithms,  which 
produced  nearly  identical  time-accurate  results. 

Because  accurate  resolution  of  the  shock  in  this  time-accurate  problem 
is  of  importance,  substantial  work  was  invested  in  investigating  the  effect  of 
different  limiters  in  resolving  the  shock  structure.  Figure  7  shows  the  density 
near  the  shock  in  the  tube,  with  Van  Albeda  and  Min-Mod  limiters,  the  last 
of  which  gave  sharper  resolution.  Similar  comparisons  were  done  for  both 
the  expansion  wave  front  and  contact  surface,  shown  in  Figures  8  and  9. 

The  GASP  time-accurate  solution  was  also  compared  to  a  completely 
independent  Flux-Corrected  Transport  (FCT)  solution,  at  an  elapsed  time  of 
0.951  ms  (where  the  FCT  data  was  available).  Figures  10-13  show  the 
comparison  between  density,  velocity,  temperature,  and  pressure  for  the  two 
independent  codes,  respectively.  Once  again,  agreement  was  excellent. 
Finally,  for  comparison,  a  viscous  solution  was  run  in  the  shock  tube;  it 
compared  favorably  with  the  inviscid  results  and  demonstrated  suitable  time- 
accuracy. 

To  validate  time  accuracy  for  external  flows,  GASP  was  used  to 
calculate  properties  over  the  classic  24-degree  wedge  solution  of  Holden  at 
Mach  14.1  and  Reynolds  number  2.36  x  lO^.  This  geometry  and  the  associated 
data  has  been  used  in  the  past  to  validate  many  other  codes.^,  and  contains 
features  which  will  be  characteristic  of  a  hypersonic  inlet  ramp.  A  grid  was 
generated  with  121  x  61  points  over  a  0.439  m-long  wedge  surface,  as  is 
depicted  in  Figure  14. 


^Rudy,  D.,  Thomas,  L.,  Kumar,  A.,  Gnoffo,  P.,  and  Chakravarthy,  S.,  "A  Validation  Study  of 
Four  Navier-Stokes  Codes  for  High  Speed  Flows,"  AIAA  89-1838,  AIAA  20th  Fluid  DYnainics, 
Plasmadynamics,  and  Lasers  Conference  Buffalo,  N.Y.  June  1989. 
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Figure  14 


Figures  15,  16,  and  17  show  density,  pressure,  and  Mach  number 
contours  at  times  1  ms,  3  ms,  and  6  ms,  respectively.  Figures  18  and  19  show 
velocity  vector  plots  of  the  wedge  flowfield  at  times  1  ms,  2  ms,  and  5  and  6 
ms.  The  separated  region  is  captured,  and  grows  through  the  calculation,  in 
exact  agreement  with  the  experimental  results.  Figure  20  shows  the  time 
history  of  the  surface  pressure,  which  compares  very  well  with  the 
experimental  data.  This  simple  case  is  obviously  far  removed  from  the  fully 
three-dimensional  inlet,  yet  it  serves  to  provide  confidence  in  the  code’s 
ability  to  handle  that  flow  with  time  accuracy. 

Geometry  Selection: 

An  important  element  of  this  work  was  the  selection  of  a  suitable  generic 
scramjet  inlet  geometry  with  realistic,  three-dimensional  effects,  and  with  the 
potential  for  application  to  a  realistic  engine  system.  This  work  is  not  tied  to 
any  particular  engine  design  or  vehicle  concept  per  se,  so  it  was  desired  to  find 
a  geometry  which  would  exhibit  physical  results  of  interest,  while  at  the  same 
time  being  as  simple  as  possible  to  generate  a  computational  grid  structure. 

For  this  work,  two  baseline  geometries  were  selected  Although  this 
project  was  directed  towards  a  fully-three-dimensional  solution,  initial 
studies  and  validation  were  done  in  two  dimensions.  Subsequent  studies 
used  the  fully  three-dimensional  geometry. 

Substantial  efforts  were  devoted  to  selecting  the  best  way  to  represent  a 
three-dimensional  inlet  in  two  dimensions.  The  simplest  solution  would  be 
to  extract  inlet  geometries  which  are  two-dimensional  slices  from  the  three- 
dimensional  form.  However,  that  does  not  provide  comparable 
thermodynamic  or  velocity  conditions.  For  instance,  if  the  three-dimensional 
inlet  has  sidewall  compression,  a  two-dimensional  slice  will  have  less 
comression,  and  therefore  not  reproduce  a  two-dimensional  slice  of  the  three- 
dimensional  flow.  As  a  result,  numerical  tests  will  be  done  on  two  separate 
sets  of  inlets,  one  for  two-dimensions  and  another  for  three.  Initial  work  has 
been  done  on  a  two-dimensional  configuration,  for  code  development  and 
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Figure  20 


validation  purposes.  Subsequent  studies  will  be  conducted  on  three- 
dimensional  geometries. 

For  the  two-dimensional  work,  two  seperate  geometries  were  selected. 
The  first  is  the  Rensselaer  Polytechnic  Institute  generic  scramjet  model.  This 
model  has  been  tested  by  Minucci  and  Nagamatsu,  across  the  Mach  number 
range  to  18,  and  so  data  is  available  for  comparison^.  Another  geometry  was 
designed  in-house  to  provide  shock-on-lip  flow  conditions  at  design,  to 
minimize  shock-boundary  layer  interactions  inside  the  combustor  channel. 
This  shape  was  adopted  because  the  RPI  model  does  not  provide  such 
operating  conditions,  and  the  computational  solutions  were  exhibiting 
extensive  separated  regions  on  the  combustor  walls  which  were  deemed 
unrealistic  for  practical  engine  designs. 


Since  the  ultimate  goal  of  this  program  was  to  study  unstart  physics 
with  a  fully-three-dimensional  geometry,  the  NASA  Langley  "Parametric 
Scramjet  Module"  was  selected  for  the  bulk  of  these  studies,  as  pictured  in 
Figure  21.  This  design  has  several  advantages,  primarily  in  that  there  is 
unclassified  experimental  data  available  for  code  validation  (though  not 
including  unstarts),  and  it  is  fully  three-dimensional,  including  the  use  of 
sidewall  compression  ramps.^  This  design  also  includes  injection  struts  in 
the  main  flow.  It  was  hoped  that  such  struts  could  eventually  be  included  in 
the  computational  grid,  but  that  was  not  achieved  during  the  reporting 
period. 


Grid  Generation 


2  Minucci,  M.,  and  Nagamatsu,  H.,  "Experimental  Investigation  of  a  2-D  Scramjet  Inlet  at  Mach 
Numbers  of  8  to  18  and  Stagnation  Temperatures  of  4,1(X)K,"  AIAA  91-0013,  29th  Aerospace 
ScierKes  Meeting,  Reno,  Nev.  January  1990. 

^Trexler,  C.,  and  Souders,  S.,  "Design  and  Performance  at  a  Local  Mach  Number  of  6  of  an  Inlet 
for  an  Integrated  Scramjet  Concept,"  NASA  Tedmical  Note  TN  D-7944,  August  1975. 
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Figure  26  presents  the  velocity  vectors  and  pressure  distribution  in  the 
separated  region  at  steady  state.  All  of  the  features  of  the  shock/boundary 
layer  interaction  have  been  captured  correctly.  Figure  27  shows  the  time 
history  of  the  evolution  of  the  separation  region;  note  the  rapid  growth  of 
this  zone.  Corresponding  velocity  vectors  are  shown  in  Figure  28,  with 
snapshots  at  1,  1.5,  and  2  ms  after  initial  shock  impingement.  The 
corresponding  density  contours  are  shown  in  the  following  Figure  29. 


Computational  Results 

Activities  for  the  second  year  were  focussed  on  calculations  in  the 
chosen  3-D  grid,  as  shown  in  Figure  30.  This  is  a  41x61x91  3-D  grid,  with  30- 
degrees  leading  edge  sweepback.  Calculations  began  with  the  intention  of 
establishing  steady-state  solutions  at  moderate  mach  number.  Figure  31 
presents  Mach  contours  calculated  with  an  Euler  solution  at  Mach  6.  The 
fully-three-dimensional  shock  structure  has  been  captured  in  this  calculation. 

Figure  32  presents  density  contours  at  the  first  and  last  z-axis  stations  in 
the  scramjet.  Temperature  contours  are  shown  in  Figure  33,  also  at  the  first 
and  last  z-axis  positions.  Normalized  pressure  and  density  are  shon  in  Figures 
33  and  34,  where  the  effect  of  each  shockwave  in  the  engine  is  apparent.  One 
point  of  note  is  the  rapid  jump  in  pressure  and  density  corresponding  to  the 
second  shock  reflection  in  the  engine  channel,  likely  source  of  an  adverse 
shock-boundary  layer  interaction  in  a  fully-viscous  calculation. 

The  three-dimensional  nature  of  the  steady-state  Euler  solution  is 
demonstrated  by  Figure  36,  which  rpesents  the  density  distribution  at  the  z- 
axis  extreme  from  the  previous  figures.  Note  the  coalescence  of  the  effects  of 
reflected  shocks  into  a  single  adverse  pressure  rise  at  x=0.14,  located  just 
behind  the  throat  of  this  chosen  geometry. 

These  inviscid  Euler  calculations  are  contrasted  against  the  fully 
viscous  Navier-Stokes  solutions  presented  in  Figure  38.  Note  for  instance  the 
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Velocity  vector  plot  in  the  combustor  region 
at  various  elapsed  times 
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Figure  30:  Computational  Grid  for  3-D  Scramjet  Engine 
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Figure  38:  Navier-Su^es  Solution 
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density  contours  at  three  z  axis  locations,  and  the  noticable  differences 
between  these  solutions  and  the  Euler  results.  Most  important  is  the 
development  of  a  large  separation  zone  on  the  initial  inlet  surface,  eminating 
from  the  shock-boundary  layer  interaction  off  the  first  reflected  shockwave. 
This  separation  zone  was  seen  to  continue  to  grow  throughout  the  time- 
accurate  solution,  making  establishment  of  a  steady-state  result  difficult.  Note 
that  this  calculation  was  done  assuming  a  laminar  boundary  layer,  so  the 
effects  of  boundary  layer  separation  are  undoubtedly  magnified  over  those 
expected  with  a  turbulent  boundary  layer. 

Velocity  contours  presented  in  Figure  40  show  a  similar  picture  of 
separated  flow  with  laminar  flow.  Note  the  vertical  nonuniformity  of  the 
separated  region,  as  demonstrated  by  the  three  z-axis  locations  plots.Velocity 
vectors  and  streamlines  are  shown  in  the  following  figure.  Figure  42  presents 
density  contours  through  the  y-axis,  along  with  temperature,  x-axis  velocity 
component,  and  streamlines  in  Figures  43-45. 

A  refined  grid  was  used  repeat  these  calculations.  A  41x71x121  grid, 
still  at  sweep  angle  of  30-degrees,  was  selected,  as  shown  in  Figure  46.  Note 
that,  in  comparison  to  results  shown  in  Figure  36  (the  low-density  grid). 
Figure  47  presents  density  levels  of  comparable  shape  and  magnitude,  though 
details  of  shock  interactions  are  more  clearly  resolved  on  the  higher-density 
grid.  Unfortunately,  fully-time  accurate  3-D  viscous  solutions  could  not  be 
run  with  this  refined  grid  in  a  reasonable  calculation  time. 

The  next  and  remairung  phase  of  this  project  would  be  the  inclusion  of 
turbulence  in  the  calculatiorw,  to  correct  the  highly-separated  flow  region  and 
avoid  the  exer-expanding  separation  zone  observed  with  the  laminar-flow 
solutions.  Next,  with  this  demonstrated  time-accurate  calculation  capability,  a 
steady-state  starting  condition  would  be  introduced  downstream  inside  the 
engine. 


Figure  40 :  Navier-Stokes  Solution 

U-Velocity  Contours  at  various  Z  locations 


Figure  A4 :  Navier-Stokes  Solution 

U-Velocity  Contours  at  various  Y  locations 
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Figure  45 :  Navier-Stokes  Solution 

Velocity  Vectors  and  Streamlines  at  various  Y  locations 
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Figure  46:  Computational  Grid  for  3-D  Scramjet  Engine 

Leading  Edge  Sweep  =  30° 

Grid  Dimensions:  41  x  71  x  121 


Figure  54 :  Navier-Stokes  Solution 

U-Velocity  Contours  at  various  Y  locati(Mis 
(41x71x121) 


